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Abstract -We show that low density homogeneous phases of self propelled hard disks exhibit a 
transition from isotropic to polar collective motion, albeit of a qualitatively distinct class from 
the Vicsek one. In the absence of noise, an abrupt discontinuous transition takes place between 
the isotropic phase and a fully polar absorbing state. Increasing the noise, the transition becomes 
continuous at a tri-critical point. We explain all our numerical findings in the framework of 
Boltzmann theory, on the basis of the binary scattering properties. We show that the qualitative 
differences observed between the present and the Vicsek model at the level of their phase behavior, 
take their origin in the complete opposite physics taking place during scattering events. We argue 
that such differences will generically hold for systems of self-propelled particles with repulsive 
short range interactions. 


S 


Self-propelled particles borrow energy from their envi- 
gronment and convert it to translational motion. Depend- 
Qing on their density and interaction details, energy taken 

I-bp on the microscopic scale can be converted into macro- 

^_^scopic, collective motion. A transition from an isotropic 
^ phase to a polar, collective phase is observed. Most theo- 
(N retical knowledge about such active matter and its transi- 
tions has been developed for Vicsek-like models: point par- 
tides travel at constant speed, to represent self-propulsion. 


and their interaction rules comprise both explicit align¬ 
ment and noise, to account for external or internal pertur¬ 
bations. In addition, self-diffusion noise may be incorpo- 
rated in the individual motion. Typically, the alignment 
strength and noise intensity control the transition between 
I!* the isotropic and the polar phases 


(N 


• ^ In order to describe more closely experiments with bac- 
teria [7||^ or with self-propelled actin filaments , other 
^ models incorporates the alignment mechanism through the 
steric interaction of elongated objects [TT|jT^ . The inter¬ 
action accounts for the body-fixed axis of these particles, 
which is in some sense more microscopic than the purely 
effective (mesoscopic) alignment rules in Vicsek-like mod¬ 
els. Steric interaction is not the only microscopic interac¬ 
tion that can lead to collective behaviour. An alternative 
way to use dynamics for aligning (hard or soft) circular 
disks is to incorporate inelasticity or softness in the inter¬ 
action rules [l4}jr^ . Not surprisingly, also these models 
exhibit a transition to a globally aligned phase. 


A model experimental system of self-propelled disks has 
been realized by vertically vibrating millimeter-sized circu¬ 
lar metal disks with two different “feet” (uni- When iso¬ 
lated, such particles advance in the direction of their body 
axis, with a reasonably high persistence length. When 
sufficiently crowded, they start to align to each other and 
move along in dense clusters. Successful modelling of this 
experiment, in the form of Newton’s equations for rigid 
disks with a force acting on the center and a torque turn¬ 
ing it around, plus slightly inelastic collisions and some 
active noise included, have demonstrated that purely re¬ 
pulsive mechanical interaction with no explicit alignment, 
were sufficient ingredients to observe a transition to col¬ 
lective motion [^. The nature of this effective alignment 
remains however mysterious. Is inelasticity a required in¬ 
gredient? Or, as suggested in [T8||^ , does it come from 
repeated recollisions of particle pairs, during which the ve¬ 
locity vector repeatedly converges against the body axis, 
but also the body axis has time to relax towards the time- 
average of the velocity axis. Given that the experiment is 
very much the one closest to a simple hard-disk liquid and 
lends itself to a description in terms of statistical physics, 
understanding precisely its phase transitions, the effective 
alignment mechanism, and the role of the noise would be 
a major step towards a theoretical description of active 
liquids in general. 

In this letter we provide such an understanding, re¬ 
stricting ourselves to the analysis of spatially homogeneous 
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Figure 1: (a) An isolated self-propelled particle converges to 
its stationary state where velocity v and polarity h are parallel, 
(b) A single binary “scattering event” can consist of many hard¬ 
disk elastic collisions, (c) Stable phases in the absence of noise. 
Between isotropic ((V^) ~ 0) and polar (('0) = 1) phase is a 
discontinuous transition. 


phases in the low density limit. To do so we proceed in 
three steps, (i) We perform molecular dynamics of the 
model equations for purely elastic interactions, with and 
without noise: In the absence of noise, the system exhibits 
a strongly first order transition from the isotropic to the 
collective motion phase (see Fig. HD- Above a finite level 
of noise, the transition becomes second order - a tricritical 
point exists. This establishes the phase behaviour which 
we will explain from theoretical considerations, (ii) We 
analyze the model equations on the grounds of the Boltz¬ 
mann equation, by making use of a recently proposed ob¬ 
servable (p • ^p) which quantifies the non-conservation of 
momentum . This observable allows to span the bridge 
from the microscopic dynamics, in particular binary col¬ 
lisions such as depicted in Fig. [^, to macroscopic order 
parameters. From a direct sampling of all possible binary 
scattering events, we obtain an excellent quantitative pre¬ 
diction of our numerical findings, (iii) We scrutinize the 
very peculiar dynamics of a collision between two self pro¬ 
pelled disks and explain the specific shape of the scattering 
function that was obtained numerically in (ii). We further 
find that recollisions are not necessary for the observed 
alignment, contrary to our previous belief. 

Model of self-propelled hard disks. The model consists 
of N hard disks in a square box of size Lx with periodic 
boundary conditions. The density is p = Njl?. Particles, 
being self-propelled, relax to a stationary speed uq. As 
units of length and time we choose the diameter do of the 
particles and do/uo, respectively. A particle i has coordi¬ 
nates r^, velocity v^, and a body axis given by the unit 
vector fi^ (see Fig. Ef)- Between collisions, it evolves ac¬ 
cording to the equations 

= Vi, (la) 

Tv^-Vi = ni - Vi , (lb) 

= {Hi X Vi) X Hi . (Ic) 


The competition between the self-propulsion n and the vis¬ 
cous damping — v in Eq. (lb) lets the velocity relax to n 
on a timescale r^. Similarly, in Eq. (Ic), the polarity n un¬ 
dergoes an overdamped torque that orients it toward v on 


a timescale r^. Interactions between particles are elastic 
hard-disk collisions which change v but not n. After such 
a collision, v and n are not collinear, and the particles un¬ 
dergo curved trajectories which are either interrupted by 
another collision (Eig. [If)> or the particles reach their sta¬ 
tionary state, where v = n and the trajectory is straight 
at a speed uq = 1 (Eig. ED- The final direction of v (equal 
to that of n) depends on the parameter 

a = (2) 

which can be understood as the persistence of the po¬ 
larity n. Linearizing the evolution equations around the 
stationary state, one can show that the final polar an¬ 
gle is given by the weighted average of the initial angles, 
{On + a0y)/{l + (a). When a ^ 1, n is practically always 
directed along v. 

On top of the deterministic trajectories given by the 
Eqs. 0, we add some angular noise by the following pro¬ 
cedure. Given a time step St 0, we rotate and Ui 
by the same angle r]i{t), distributed normally with zero 
mean and variance 2DSt^ where the constant D > 0 fixes 
the level of the angular noise. Noises of different particles 
are statistically independent. We choose St much smaller 
than all other timescales in the dynamics. The relevant 
parameter to characterize the angular noise is then D/X, 
where A = 4p/7r is the characteristic scattering rate of the 
system, which is proportional to the density [^ . 

Moleeular dynamies (MD) simulations. We now es¬ 
tablish the phase behaviour of the model for N particles. 
MD simulations were performed at = 4 with N = 1000 
or N = 4000, focusing on the dilute regime p 1 (see 
below for a discussion of the effect of Ty). We are thus 
left with two microscopic parameters, namely a and D/X. 
Also, the system size is chosen not too large, in order to 
keep the system spatially homogeneous, which we have 
checked by visual inspection. We measured the order pa¬ 
rameter f;{t) = \"^i^i{t)\/N, which is of order l/\fN for 
the isotropic state and of order one for the polar state. 

Let us first look at the case without angular noise, 
D/X = 0. We initialized simulations from random iso¬ 
tropic conditions and waited for the isotropic state to even¬ 
tually destabilize. When a stationary state was reached, 
we started to average the order parameter over time, {fj). 
As shown in Eig. [^, we found the isotropic state to be 
stable at low values of a, whereas it becomes unstable 
at larger values, in favour of a polar state. Between the 
two phases, an abrupt discontinuous transition takes place 
at . Quite remarkably, in the whole polar phase the dy¬ 
namics converges to = 1, where particles are all strictly 
parallel. Eurther, choosing some random state with ^ 1 
as initial condition, we found that the polar state '0 = 1 
is stable for all a > 0, in particular also when a < a^. 

In Eig. [^, we show again the (now rescaled) order pa¬ 
rameter in the isotropic state, this time for different densi¬ 
ties. Eor a given density, the data for different values of N 
collapse, showing that finite-size effects are under control. 
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Figure 2: (a) Order parameter in the isotropic phase, without 
angular noise. N = 1000 (circles) and 4000 (crosses). From 
right to left: theory at p ^ 0, p = 10“^, 0.01, 0.02, 0.05, 
0.1, 0.2, 0.3. (b) Dependence of the isotropic-polar transi¬ 

tion on the density, with no angular noise, (c) Transition 
lines in the (a, D/A)-plane, at fixed density p = 0.01. Solid 
lines are theoretical results at = oo. Upward triangles (red) 
and downward triangles (blue) are transitions measured nu¬ 
merically by respectively increasing and decreasing D/X quasi- 
statically. (d) Order parameter obtained by increasing (red) 
and decreasing (blue) the angular noise. 


In all cases we observed convergence to the fully polar state 
beyond the points shown. The theoretical framework used 
below adds the line for p ^ 0, drawn in black. Increasing 
density here clearly favours the polar state. The departure 
of the transition due to density effects, 1 — (a*(p)/(a*(0), is 
plotted in Fig. E>- 

Adding angular noise to the trajectories quite changes 
the picture. During simulations we first increased D/X 
quasi-statically and then decreased it again. The transi¬ 
tion was measured by looking at the maximum of the fluc¬ 
tuations of the order parameter among many realizations 
of the dynamics. The resulting phase diagram at density 
p = 10“^ is shown in Fig. [^. In agreement with intu¬ 
ition,the isotropic state is always stable at strong enough 
angular noise. When decreasing the noise we pass into the 
polar phase, but the nature of this transition can be either 
discontinuous or continuous. For values a < ac ~ 0.157 
the transition has some hysteresis, as indicated in the up¬ 
per panel of Fig.[^. The discontinuous nature of the tran¬ 
sition is thus robust when adding angular noise, and the 
phase areas in Fig. overlap, presenting an area which 
we could call coexistence region if the system were not 
homogeneous. For a > o^c, the hysteresis is no longer 
observed at our level of numerical precisiorQ as can be 
seen in the lower panel of Fig. ii. At Q^c, the coexistence 

^In the last three points of Fig. [^, for a > 10 there is a slight 
hysteresis which disappears for slower annealing rate of D/X. 


zone vanishes into a single line of transition (tricritical 
point p2||^ ). Interestingly, once in the polar phase, fur¬ 
ther increasing a leads to a re-entrant transition towards 
the isotropic phase. 


Kinetic theory framework. We now rationalize these 
numerical observations in the context of kinetic theory, 
using the properties of the binary scattering, following 
Refs. nai . We summarize here the procedure described 
in Ref. IT In the dilute regime, the mean free-flight time 
is long enough so that particles have mostly reached 
their stationary velocity vq before interacting with another 
particle {ry^Tn A“^). The binary scattering of self- 
propelled particles does not conserve momentum, and that 
is why a polar state can emerge from an isotropic initial 
condition. Assuming molecular chaos and choosing the 
von Mises distribution as an ansatz for the angular distri¬ 
bution of the velocities, one can write down an evolution 
equation for the order parameter f;. This equation can 
then be expanded, up to order to study the stability of 
the isotropic phase 


- (m- 

where the coefficients are given by 


M := (p-'5p)o> 

^ := (d -cosA)p-(5p)g, 
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In the stationary state, the left-hand size of Eq. (§ van¬ 
ishes. The transition line is obtained by solving the equa¬ 
tion p{a^) = D/X for a*, while the sign of f{a^) at the 
transition tells whether it is continuous or discontinuous. 
The coefficient p is exact within the assumptions of kinetic 
theory, while f should depend on the ansatz used for the 
angular distribution. Both coefficients are an average over 
all pre-scattering parameters, as given in Eq. (|^, where 
b is the impact parameter and A is the angle between the 
incoming particles’ velocities. The averaging needs not be 
done over the norms of the velocities, since those are fixed 
to 'T’o = 1- We have checked explicitly that this assumption 
holds very well in the numerical simulations in the dilute 
regime. In Eqs. 0 and 0, p is the pre-scattering mo¬ 
mentum of the two colliding particles, and 6 p is the change 
of their momentum by the scattering event. In contrast 
with Ref. [^, the above equations explicitly identify the 
forward component of the momentum change p • (ip, as 
the proper quantity to describe the alignment in a binary 
scattering event. The predictions thus depend only on the 
microscopic details of the deterministic dynamics through 
the scattering function p • Sp{b^ A), which we now explore 
for the model 0- 

Binary scattering. Remind that a scattering event can 
consist of several, if not many, hard-disk collisions, such 
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Figure 3: Top: Colour maps of the full scattering function p • (5p(5, A). The geometry of the collisions and the colour scale are 
shown in the bottom-left panel. Bottom—middle: the partially integrated scattering function as a function of the incoming 
angle A, for different values a, = 4. Right: Fully integrated scattering functions fi and defined in Eqs. @ and 0 plotted 
as a function of a for different Ty. The red vertical line indicates the transition in the absence of noise. 


as depicted in Fig.[^. Before the first collision, both par¬ 
ticles are taken to be in the stationary state, i.e. v = h. 
After the collision, we integrate Eqs. 0 numerically un¬ 
til another collision possibly occurs. The binary scat¬ 
tering is considered over when both particles have again 
reached their stationary state and are heading away from 
each other. Repeating the procedure in the whole range 
of initial parameters (6, A) yields the scattering function 
p • Sp{b^ A). Figure shows this function for several a, 
as well as its integral over b and the full integrals yielding 
the coefficients /i and 

Let us stress that p • ^p is not changed by the colli¬ 
sion itself, which conserves momentum. All (dis) alignment 
must here come from the relaxation of the post-collisional 
value of |v| to unity. Fig.[^-d shows that scattering at low 
angular separation, small A, always creates forward mo¬ 
mentum. In other words, two nearly parallel particles that 
interact become even more parallel, which gives rise to an 
effective alignment, p • ^p > 0. On the other hand, for 
small enough a, particles that enter in interaction frontally 
(A ^ tt) tend to disalign, except for special symmetry such 
diS b ^ 0. Increasing a favours aligning scattering events 
until eventually only aligning events remain. This is best 
summarized by integrating out all parameter dependence 
except the incoming angle, as is plotted in Fig. Il,f. 

The coefficients /i and ^ are then obtained by integrat¬ 
ing over A, with weights prescribed by Eqs. 0 and 0. 
Their dependences on the microscopic parameters of the 
dynamics, a and r^, are shown in Fig. [^,h. In the ab¬ 
sence of noise, the transition occurs for a = such that 
= 0; is negative, hence the discontinuous tran¬ 

sition. When angular noise is added, the transition is ob¬ 
tained by solving the equation /i((a*) = D/X. From the 
shape of the curve one obtains two values a±, with 

a- and ^ oo when D ^ 0. As D/A is increased, 

^{a-) eventually becomes positive and the transition turns 


continuous at a tricritical point {ac^Dc/X). Note that 
^((a+) > 0: the re-entrant transition is always continuous. 
Regarding the role of r^, a- is practically independant on 
r^, while increases when decreases. The theoretical 
predictions are shown as solid lines in Fig. 0- The agree¬ 
ment with the MD simulations data for density p = 10“^ 
is excellent. The small shift of the measured transition 
lines to the left with respect to the theoretical one comes 
from finite-density effects. 

Finally, we also learn from the examination of the scat¬ 
tering maps that, in the absence of noise, the polar phase 
^ = 1 is actually an absorbing phase [^: this is because 
all binary scattering events at small A have p • (^p > 0. 
When all particles in a system are sufficiently parallel, bi¬ 
nary scattering events can only align the system more. 
This is true for all a, and most remarkably for a ^ 0. 

Altogether, our kinetic theory description, using the von 
Mises ansatz for the angular distribution, captures quanti¬ 
tatively all the phenomenology reported in the numerical 
simulations at low enough density. It however relies on the 
numerical evaluation of the scattering maps. In the last 
part of the paper, we would like to provide some intuition 
on the origin of the peculiar form of these maps. Also, we 
will elucidate the role of the multiple collisions which can 
take place during a scattering event. 

Analytical limits. When we take the limit ^ 0, 
the vector n has no persistence at all. After two parti¬ 
cles collide elastically, the n rotate to their respective v 
instantaneously, the two particles cannot collide a second 
time and the post-collisional relaxation of |v| to unity oc¬ 
curs on a straight trajectory. We take advantage of this 
simplification to compute the scattering map in the case 
a = 0. Because velocities have equal modulus, we can use 
the nice visualisation of an elastic collision as a rotation 
in the center-of-mass frame (Fig.[^,b). One can then find 
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Figure 4: (a,b) Elastic collisions are rotations by an an¬ 

gle 2(3 — 2 asin(6) in the reference frame of the center-of-mass 
velocity v. Primed quantities (green) are post-collisional, black 
ones are pre-collisional. (c) The scattering function p • (5p re¬ 
sulting from a single collision, (d) A cut through (c) for the 
value of b used in (a,b). 


an analytic expressiorj^for p • Sp{b^ A), which is plotted in 
Fig#. This plot nicely completes the series of varying a 
from Fig. [^-d. The essential structure of the scattering 
maps can be accounted for by the sole ingredients present 
in the case a = 0. Inelastic collisions and persistence of n 
are not essential. Conversely, non-conservation of momen¬ 
tum due to the relaxation of the particle speeds to unity 
is crucial. 

At values a > 0, the vector n has non-vanishing per¬ 
sistence, which results in curved trajectories and possi¬ 
ble recollisions. This could lead to the effective align¬ 
ment mechanism which was proposed to be at the root 
of the collective motion in the experiment of vibrated po¬ 
lar disks [^- 20 . Concerning the influence of recollisions 
in the scattering, we ask whether they contribute signifi¬ 
cantly to fi and In particular, how numerous are they 
depending on the scattering geometry (6, A), what is their 
statistical weight and how much do they affect the scat¬ 
tering map p • Sp{b, A)? One can get an intuitive picture 
of the recollision mechanism looking at a simple case, the 
symmetric {b = 0) binary scatterings in Fig. |^. A first 
collision makes v rotate from being equal to n to point¬ 
ing away from the other particle. On the following trajec¬ 
tory, V and n relax towards each other, on their respective 
timescales Ty and r^. For small enough Ty, the persistence 
of n allows for a recollision to take place, and both vec¬ 
tors get closer to each other from one collision to the next. 
The highly symmetric case depicted in Fig. ii is however 
so degenerate that it has either one or infinite collision^ 


^p • Sp{b, A) = 2 cos ^ 


E 


-cosf 


^ l+2sb\/l-b^ sin A 

^In the symmetric case, linearizing the dynamics around the sta¬ 
tionary point provides analytical results about this mapping from 



Figure 5: Recollisions, (a) Symmetric binary scatterings for 
0 = 5. (b, c) Number of recollisions in each scattering (top 
rows) and their relative contribution to p • (5p (bottom rows), 
o = 0.1 for (b), o = 10 for (c). 


Unfortunately, the analytic treatment of asymmetric col¬ 
lisions is much too technical to gain physical insight from 
it and we prefer coming back to the numerical sampling 
of the scattering events. 

For a = a_, as shown in Fig. HI (top row), recolli¬ 
sions can be numerous, but take place only in a small re¬ 
gion of the parameter space (6, A), the smaller the larger 
Ty . We measure their relative contribution to p • (^p by 
T = i sin(A/2)(p • (5p - p • (5p|i)/|p • (5p|, where p • (5p|i 
is the contribution of the first collision only. As shown 
in Fig. (bottom row), even for the smallest value 
Ty = 0.04, this contribution is locally less than 5 %. For 
larger r^, the contribution is vanishingly smaller. Thus, 
recollisions has no practical influence on the transition at 
a = a-. Conversely, as shown in Fig. (top row), for 
a = the number of recollisions is not so large, but they 
occur for wider ranges of scattering parameters and their 
relative contribution to p • (^p is significant, see Fig. 
(bottom row). This is all the more true for smaller Ty and 
is responsible for the dependance of fi and ^ on Ty in this 
regime. 

Discussion. Let us recast our main findings and the 
understanding we obtained from the above analysis. A ho¬ 
mogeneous system of self propelled hard disks, obeying 
the minimal deterministic Eqs. exhibits a sharp dis¬ 
continuous transition from an isotropic to an absorbing 
polar collective motion state, when increasing the persis¬ 
tence of the body axis vector n. Adding angular noise, the 
transition becomes continuous via a tricritical point. At 

one collision to the next. For o: < 1, there is no recollision. For 
a > 1, an infinite number of recollisions occur, with a collision rate 
that converges to a positive constant {a < 2) or diverges (a > 2). 
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fixed noise level, increasing further the persistence of the 
body axis vector, a re-entrant continuous transition from 
the polar state to the isotropic state occurs. 

All these macroscopic behaviours can be understood by 
the study of the scattering maps p • which describe the 
level of alignment of all possible scattering events. The 
most striking feature of these maps is that scattering at 
low angular separation always creates forward momentum. 
This is even true for a = 0, for which, we could prove 
it explicitly from a mechanical analysis of the collisions. 
This feature explains the stability of the polar phase in the 
absence of noise. For small a, the dis-alignment, present 
in the frontal collisions is strong enough to stabilize the 
isotropic state. For a > this does not hold anymore, 
and the isotropic state becomes unstable. The coexistence 
of stability of the two states for a < guarantees the 
transition to be discontinuous. 

The effect of noise is most easily understood from the 
shape of the curves ii{a) and ^(of). The key point here 
is that in the present scheme of approximation - Boltz¬ 
mann equation plus von Mises Ansatz - the computation 
of ^ is identical to that of /r, apart from the additional 
asymmetric factor (1/2 — cos A). This does not rely on 
the specificity of the microscopic model. As a result, the 
curve is similar to the curve with a shift towards 
larger a. Adding noise simply shifts the transition towards 
larger a\ at some point ^ changes signs at the transition, 
which becomes continuous, hence the tricritical point. 

Finally, the re-entrant transition from polar back to iso¬ 
tropic at large a (while increasing a) is a direct conse¬ 
quence of the non-monotonous shape of At moder¬ 

ate (a, the effective alignment /i increases with the persis¬ 
tence of n. However, a too large persistence reduces the 
alignment, because the scattering event does not reori¬ 
ent n significantly, and no memory of the collision is kept: 
V relaxes to practically the same n as before the scatter¬ 
ing. Some alignment still takes place but is most easily 
destroyed by small amounts of noise. Interestingly, while 
the recollisions play no role at the first transition from the 
isotropic to the polar state, they here increase the effective 
alignement and thereby postpone the re-entrant transition 
to larger a. 

As can be seen from the above discussion, the results 
do not depend on the details of the scattering maps, as 
long as binary scattering at low angles produces effective 
alignment. The partially integrated scattering function 
(Fig. |^,f) is then positive for small angles. Models which 
fall into this class are (i) the here described hard disks 
with V and n, (ii) the inelastically colliding hard disks with 
only V [^, and (hi) the soft disks described in Ref. pT) . 
All three are self-propelled in the sense that their free- 
flight speed relaxes to uq. 

It is important to note that the distinctive behaviour of 
the above class of system is markedly different from the 
original Vicsek model, in which binary scattering at low 
angle effectively disaligns. In the Vicsek model, the par¬ 
tially integrated scattering function Fig. 2a] behaves 


precisely in the opposite way of the one here. It is there¬ 
fore unlikely that a coarse-graining of self-propelled hard 
disks (or any of the above three models) would yield the 
Vicsek model. Altogether, the present new class of mod¬ 
els, including the self-propelled hard particles, could well 
play for the theory of “simple active liquids” the role that 
hard spheres play for the statistical mechanics of gases. 
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